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Abstract 

We have updated the parton and hadron cascade model PACIAE 2.0 (cf. 
Comput. Phys. Comm. 183 (2012) 333) to the new issue of PACIAE 2.1. 
The PACIAE model is based on PYTHIA. In the PYTHIA model, once the 
generated particle/parton transverse momentum pt is randomly sampled, 
the p x and p y components are originally put on the circle with radius px 
randomly. Now it is put on the circumference of ellipse with half major 
and minor axes of pr(l + 8 P ) and pr(l — S p ), respectively, in order to better 
investigate the final state transverse momentum anisotropy. 
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Nature of problem: PACIAE is based on PYTHIA. In the PYTHIA model, once 
the generated particle /parton transverse momentum (pr) is randomly sampled, the 
p x and p y components are randomly placed on the circle with radius of px ■ This 
strongly cancels the final state transverse momentum asymmetry developed from 
the initial spatial asymmetry. 

Solution method: The p x and p y component is now randomly placed on the cir- 
cumference of ellipse with half major axis ofpr(l + 5 P ) and the half minor axis of 
pr(l — S p ) instead of circle. 

Reasons for the new version: PACIAE is based on PYTHIA, where once the gen- 
erated particle /parton transverse momentum (px ) is randomly sampled, the p x and 
p y components are randomly placed on the circle with radius of px- This is not 
only strongly canceling the final state transverse momentum asymmetry developed 
from the initial state spatial asymmetry but also inconsistent with the ATLAS ob- 
servation of the final state charged particle transverse sphericity is less than unity 

fly- 

Summary of revisions: The main revision is executed by randomly placing p x and 

Py components of the generated particle /parton transverse momentum pt on the 

circumference of ellipse with half major axis of + S p ) and the half minor axis 

of pt{1 — 5 p ) instead of circle. 

Restrictions: Depend on the problem studied. 

Unusual features: 

Additional comments: Email addresses: zhoudm@phy.ccnu.edu.cn (D.-M. Zhou), 
yanyl@ciae.ac.cn (Y.-L. Yan). 
Running time: 

• Using the attached input file of usux.dat (where the string fragmentation is 
selected and the elastic parton-parton interactions is considered only, the 
same later) to run 1000 events for the ^fs=200 GeV Non Single Diffractive 
pp collision by 21 a. tar. gz spends 0.5 minute. 
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• Using the attached input file of usu.dat to run 10 events for the 10-40% 
most central Au+Au collisions at y/swN =200 GeV by 21b. tar. gz spends 5 
minutes. 

• Using the attached input file of usu.dat to run 10 events for the 10-40% 
most central Au+Au collisions at y/sNN=200 GeV by 21c. tar. gz spends 17 
minutes. 

1. 

The large azimuthal anisotropy (the large second harmonic coefficient 
v 2 ) of the emitted particle is an important feature of the hot and dense 
medium created in the ultra-relativistic nuclear collisions. This large v 2 has 
contributed to the observation of a strongly coupled quark-gluon plasma 
(sQGP) in the nucleus-nucleus collisions at the RHIC energies 0, S 0, 0] • 

The nuclear overlap zone created in a nucleus-nucleus collisions at a given 
impact parameter possesses an almond-like spatial asymmetry. Because of 
the strong parton rescattering, the local thermal equilibrium and asymmetric 
pressure gradient may build up in this initial fireball. The asymmetric pres- 
sure gradient then drives a collective anisotropic expansion. The expansion 
along the almond minor axis (along the large pressure gradient) is faster than 
the one along the major axis. This results in a strong asymmetric transverse 
momentum azimuthal distribution and hence a large elliptic flow coefficient 
v 2 of the final hadronic state. 

As mentioned in [6], PACIAE is a parton and hadron cascade model for 
the ultra-relativistic nuclear collisions and is based on PYTHIA |7|. In the 
PACIAE model, a nucleus-nucleus collision is decomposed into a sequence 
of nucleon-nucleon (NN) collisions according to the collision geometry and 
the NN total cross section. Each NN collision is performed, in turn, by 
the PYTHIA model with the string fragmentation switched-off temporarily 
and the diquark (anti-diquark) broken into quark pairs (anti-quark pairs) 
randomly. The parton rescattering then proceeds. This parton evolution 
stage is followed by the hadronization at the moment of partonic freeze-out 
(exhausting the partonic collisions). The Lund string fragmentation regime 
and/or phenomeno logical coalescence model is provided for the hadroniza- 
tion. Then the rescattering among produced hadrons is dealt with the usual 
two body collision model j6|. 

In the PYTHIA model j?} once the transverse momentum of a final 
state hadron generated from the string fragmentation and/or the unstable 
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Figure 1: (Color on line) The charged particle i>2(pt) at mid-rapidity (\rj\ < 1.) in 10-40% 
most central Au+Au collisions at y / ijv r /v=200 GeV. The solid symbols are the STAR data 
taken from [l3| and open symbols are the PACIAE results calculated with C—2, 3, and 4. 



particle decay is randomly sampled, the p x and p y components are randomly 
placed on the circle with radius of pr- This p x and p y determination method 
may completely cancel the final state transverse momentum anisotropy de- 
veloped from the initial spatial asymmetry. The charged particle transverse 
sphericity [J, @, ^ may equal to one (isotropic). This is inconsistent with 
the experimental observation that the charged particle transverse sphericity 
is less than unity [if. Therefore we are randomly placed the generated final 
state hadrons on the circumference of ellipse with half major axis of pr(l+8 p ) 
and the half minor axis of pr(l — <5p) instead of circle now. This change is 
also introduced in the particle/parton production process of hard scattering, 
multiple interactions, initial- and final-state parton showers, as well as the 
adding of remnants [7j • This change is even introduced in the deexcitation of 
energetic quark (anti-quark) when the phenomenological coalescence model 
j6[ is selected for hadronization. Of course, a new px should be recalculated 
by p x and p y after this change. Then the transverse momentum distribution 
of final state hadron may be modified. However, if the deformation param- 
eter S p is less than unity (a small perturbation) the change in transverse 
momentum distribution may be weak. 

From ideal hydrodynamic calculations [loj one knows that the integrated 
elliptic flow parameter is directly proportional to the initial spatial eccen- 
tricity of the nuclear overlap zone. Therefore, if the nuclear overlap zone is 
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Figure 2: (Color on line) The charged particle V2(rf) (0.15 < pr < 2.) in 10-40% most 
central Au+Au collisions at ^/satat=200 GeV. The solid symbols are the STAR data taken 
from [13J and open symbols are the PACIAE results calculated with C=2, 3, and 4. 



assumed to be an ellipse with major axis of b = 2/2^(1 + 8 r ) and the minor 
axis of a = 2Ra(1 — 8 r ) for a symmetry nucleus- nucleus collision with the 
nuclear radius of Ra, we may assume 

8 p = C8 r (1) 

where C is an extra model parameter. C = corresponds to the original 
case of p x and p y put on the circle randomly. 

In order to calculate 8 r we first calculate the reaction plane eccentricity 
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according to the participant nucleons spatial distributions inside the nuclear 
overlap zone p] in the PACIAE simulation. In the above equation, a 2 x = 
x 2 — x 2 (the same for a 2 ) and x 2 (x) denotes an average of x 2 (x) over 
particles in a single event. The event average reaction plane eccentricity 
reads 

M = (fe^)- (3) 



On the other hand, the geometric eccentricity [12[| of the ellipse-like nuclear 
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Figure 3: (Color on line) The charged particle transverse momentum distribution in 10- 
40% most central Au+Au collisions at y / s]v"/v=200 GeV calculated by the PACIAE model 
with C=0, 2, 3, and 4: (a) in the full 77 phase space, (b) in \rj\ < 1. 



overlap zone is 



b 2 -a 2 



~ g V V 
Letting e g = e rp , one approximately obtains 

-r; 
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(4) 



5 r ~ (5) 



The calculated charged particle v^ipr) at mid-rapidity (\rj\ < 1.) in the 
10-40% most central Au+Au collision at v /s/v r /v=200 GeV is compared with 



the corresponding STAR data [13| in Fig. [TJ In this figure the STAR data 
are denoted by solid symbols: the black circles are measured with the event 
plane method (EP), red squares with Lee- Yang zero point method (L-YZ), 
and green triangles with four particle cumulant method (4 cumulant). The 
PACIAE results are given by open symbols: the black circles calculated with 
C=4, red squares with C—3, and green triangles with C=2. One sees in this 



figure that the STAR data [13J on the charged particle V2(pr) are able to 
be reproduced by the PACIAE calculations with C=3. The C=0 PACIAE 
results are too small compared to the STAR data. 
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A similar comparison for charged particle V2(rj) (0.15 < px < 2. GeV/c) 
in the 10-40% most central Au+Au collision at v /satat=200 GeV is given in 
Fig . |2j Here one sees again that the STAR data on the charged particle V2{rj) 



131 ] are able to be reproduced by the PACIAE calculations with C=3. We 
have to mention here that in Fig. 3 of Ref. 13J the PHOBOS data 14j were 
introduced to compare with the STAR data and to complement the lack of 
the STAR data in 1.5 < |r/| < 2.5 region. Because the PHOBOS data were 
measured for the 0-40% most central Au+Au collisions at the same energy 
but in the full pr phase space, it is not suitable to compare the PHOBOS 
data with the STAR data. Therefore we do not include the PHOBOS data 



14j in Fig. [2] here. 

We give the calculated charged particle transverse momentum distribu- 
tion in 10-40% most central Au+Au collisions at v /s/v r /v=200 GeV in Fig. |3j 
Figure |3] (a) and (b) are drawn for the full and partial (\rj\ < 1.) pseudo- 
rapidity phase space, respectively. In panels (a) and (b) the solid black 
circles, open red circles, open green triangles, and the blue line are calcu- 
lated with C=4, 3, 2, and 0, respectively. We see in this figure that the 
charged particle transverse momentum distribution is really not sensitive to 
the parameter C both in the full and partial 77 phase space, provided the 
deformation parameter 5 P is less than unity (a small perturbation). 

In addition, an extra switching parameter of iparres is introduced in the 
new issue of PACIAE 2.1. We assume that the iparres=0 is for the elastic 
parton-parton rescattering only and iparres=l for otherwise. 
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